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By means of neutron powder diffraction measurements in the 1.5-300 K range, we investigated the 
effect of the polar Bi 3+ ion on the magnetic ordering of the Mn 3+ ions in BiMn3Mn40i2, the coun- 
terpart with quadruple perovskite structure of the simple perovskite BiMn03. At all temperatures, 
the data are consistent with a noncentro symmetric spacegroup Im which contrasts the centrosym- 
metric one 12 /m previously reported for the isovalent and isomorphic compound LaMn3Mn40i2. 
This difference gives evidence of a Bi 3+ -induced polarization of the lattice in BiMn3Mn40i2. At 
low temperature, the two Mn 3+ sublattices of the A' and B sites are found to order antiferromag- 
netically (AFM) in an independent fashion at 25 and 55 K, respectively, similarly to the case of 
LaMn3Mn40i2. However, both magnetic structures of BiMn3Mn40i2 radically differ from those of 
LaMn3Mn40i2. Specifically, in BiMn3Mn40i2 the moments M^/ of the A' sites form an anti-body 
AFM structure, whilst the moments Mb of the B sites result from a large and uniform modula- 
tion ±Ms,b along the fo-axis of the moments Mb, oc in the ac-plane of an _E-type structure. The 
modulation is found to be strikingly correlated with the displacements of the Mn 3+ ions induced by 
the Bi 3+ ions. Our symmetry analysis of this correlation unveils a strong magnetoelastic coupling 
between the internal strain field created by the Bi 3+ ions and the moment of the Mn 3+ ions in the B 
sites. We ascribe this phenomenon to the high symmetry of the oxygen sites and to the absence of 
oxygen defects, two characteristics of quadruple perovskites not found in simple ones, which prevent 
the release of the Bi 3+ -induced strain through distortions or disorder. This result demonstrates 
the possibility of achieving a large magnetoelectric coupling mediated by internal strain in proper 
ferroelectrics and suggests a novel concept of internal strain engineering for multiferroics design. 

PACS numbers: 75.85.+t,75.47.Lx,75.25.-j,75.80.+q 



I. INTRODUCTION 

The possibility of mutually controlling electric polar- 
ization and magnetization in multiferroic materials may 
lead to new device concepts for electronics. The challenge 
is to design materials with an effective magnetoelectric 
coupling suitable for applications. A first strategy is to 
achieve a coupling between polar and magnetic ions in 
the same compound^—. This is typically the case of Bi- 
based compounds, such as the perovskite system Bi£?03, 
where the 6s 2 lone pair of the Bi 3+ ion forms an elec- 
tric dipole, while the B site is occupied by a magnetic 



ion, such as Mn 3+ or Fe 3+ — ~— . Despite intense research, 
modest magnetoelectric couplings have been hitherto re- 
ported in these compounds; this has been explained by 
the fact that the mere coexistence of polar and mag- 
netic ions does not necessarily lead to their coupling^—. 
Indeed, in these compounds, the ferroelectric order de- 
velops independent of the magnetic one {proper ferro- 
electricity) . Recently, the focus has thus been moved to 
improper ferroelectrics, where ferroelectricity is induced 
by magnetism. Among the manganese oxides relevant 
here, notable are the RMnOjLS. and i?Mn 2 0^-ii sys- 
tems (R=Y or rare earth), where a sizable tunability of 
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the polarization by applied fields has been reporte d 12 ' 13 . 
Various magnetoelectric coupling mechanisms have been 
considered, such as the spin-current interaction in spi- 
ral magnetic structures or the exchange striction in frus- 
trated magnets 1 ^— , and the debate has attracted a great 
deal of interest. 

In the absence of a firm theoretical framework, the 
possibility of a large coupling between polar and mag- 
netic ions remains open and new hints are provided by 
recent experiments. Notable is a striking enhancement 
of magnetoelectric response in epitaxial BiFcC>3 films, 
where strain drives a rhombohedral to pseudo-tetragonal 
phase transitio n 19 ' 20 . Although the mechanism of this 
enhancement remains unclear, this result suggests that 
the magnetoelectric coupling in proper ferroelectrics is 
sensitive to strain and structural distortions 2 ! To clar- 
ify this point would be important for multifcrroics de- 
sign. This is a challenging task for the Bi£?03 system 
owing to the complex structure-property relationships, 
e.g. polymorphism— modulated structure 2 ^ and sen- 
sitivity to oxygen vacancies^—. 

In this paper, we address the above point in a sim- 
pler system. By neutron powder diffraction, we have 
studied the effect of the polar Bi 3+ ion on the magnetic 
ordering in BiMri3Mri40i2 29 ' 30 , the counterpart with 
quadruple perovskite structure of the simple perovskite 
BiMnC>3. Quadruple perovskites are described by the 
general formula AA3B4O12 — ~— and we shall restrict 
to the single- valent manganese oxides ^4Mn3Mn40i2, 
where A is a three-valent ion, such as La 3 ^— , Pr— , 
Bi^ 9 .^, and Mn 3+ occupies both A' and B sites. The 
structure (see Fig. 1) consists of the same pseudocu- 
bic network of corner-sharing BO§ octahedra character- 
istic of simple perovskites ^4-BC>3, such as the aforemen- 
tioned BiMnC>3 and TbMnOs multiferroics. Our moti- 
vation is that, in spite of this similarity, quadruple per- 
ovskites display a simpler distortion pattern and no oxy- 
gen vacancie o 33 ' 34 ' 38 . Our results show a striking corre- 
lation between the internal strain field induced by the 
polar Bi 3+ ions and the magnetic ordering of the Mn 3+ 
ions in BiMn3Mn40i2- This conclusion provides a hint 
to achieve an effective magnetoelectric coupling between 
polar and magnetic ions in perovskite-likc compounds. 



A. Structural characteristics of quadruple 
perovskites AMri3Mn40i2 

As discussed in previous paper o 33 ' 34 , the above favor- 
able characteristics arise from the following features of 
quadruple perovskites. 1. Their structure is obtained by 
doubling the cubic unit cell axis of simple ABO3 per- 
ovskites; the new unit cell therefore contains two distinct 
A and A' sites at the corners and middle of the cube 
edges, respectively. This doubling is caused by a large 
buckling of the BOq octahedra driven by the Jahn- Teller 
(JT) distortion of the A' site. As a result, the Mn-O- 
Mn bond angle is drastically reduced to ~ 135°, an un- 



usually small value as compared to the ~ 160° value of 
simple perovskites, such as LaMn03, which typically dis- 
play orthoferrite-type distortions^ 9 .. 2. The large buck- 
ling lowers the coordination number of the A' site from 
twelve to four, which prevents oxygen vacancies, as these 
would destabilize the structure. Indeed, no evidence of 
vacancies has been hitherto reported in any quadruple 
perovskit o 29 ' 33 ' 34 ' 37 ' 38 . No extra oxygen atoms can be 
hosted either within the structure owing to its high- 
densit y 35 ' 37 . 3. The simple distortions arise from the 
high symmetry of the cubic 7m3 phase, where all oxygen 
atoms occupy a unique high symmetry to., site, which 
only permits a rig id tilt of the B0 6 octahedral! The 
JT distortion of the octahedra requires a cubic 7ro3 to 
monoclinic 12 /to phase transition typically realized at 
low temperature. 

B. Single- valent Mn 3+ quadruple perovskites 
LaMn 3 Mn 4 Oi2 and BiMn 3 Mn 4 Oi 2 

In view of the above considerations, quadruple per- 
ovskites are a model system to investigate the coupling 
between Bi 3+ and Mn 3+ ions and BiMn3Mn40i2 is an 
obvious choice. It was previously reporte d 29 ' 30 that the 
nuclear structure of this compound is similar to that of 
the isovalent LaMn3Mn 4 0i2 34 ' 36 , where La 3+ instead 
of Bi 3+ occupies the A site. A structural study on sin- 
gle crystals^ 9 - shows that the main difference between the 
two structures is the center of symmetry breaking in the 
former, which is described by the Im, instead of 12 /m, 
space group. This study confirms that this breaking is 
caused by the polar properties of the Bi 3+ ion, which 




FIG. 1. (Color online) Nuclear unit cell of the quadruple per- 
ovskite structure of BiMn3Mn40i2- Note that the Mn 3+ ions 
occupy two distinct A' and B Jahn- Teller sites with square 
and octahedral coordination, respectively. 
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leads to a sizable spontaneous polarization ~ 7^C cm -2 
at room temperature. In addition, BiMn 3 Mn 4 0i2 dis- 
plays an AFM transition at Tjy = 55 K concomitant to 
an anomaly of the dielectric constant, which suggests a 
coupling between AFM order and polarization. Interest- 
ingly, LaMn3Mn40i2 exhibits a similar AFM transition 
in the same temperature range T/v=78 K but no dielec- 
tric anomal y 34 i 36 . The above results motivate the present 
study of the magnetic structure of BiMn3Mn40i2- 



II. EXPERIMENTAL METHODS 

1. High-pressure synthesis of BtMn-iMriiOu powders 

Powder samples of single-phase BiMn3Mn40i2 were 
synthesized under high pressure in a multi-anvil appara- 
tus at 4 GPa and 1000 °C for one hour as discussed in 
detail elsewhere 2 ^. After the high pressure and high tem- 
perature synthesis, the samples were quenched down to 
room temperature in less than 1 minute in order to sta- 
bilize the high pressure phase. In order to check phase 
purity, the as-prepared samples were first analyzed us- 
ing a commercial Siemens X-ray powder diffractometer 
equipped with a Cu K Q radiation. 



2. Neutron powder diffraction measurements 

The nuclear and magnetic structures were studied at 
300 K and in the 1.5-100 K range by means of neutron 
powder diffraction data collected at the 3T2 and G4- 
1 diffractometers, respectively, of the Laboratoire Leon 
Brillouin in Saclay, France. The 3T2 instrument is a 
high-resolution diffractometer with wavelength A=1.225 
A, whilst G4-1 is a high intensity diffractometer opti- 
mized for magnetic structure determination with A=2.43 
A. Both nuclear and magnetic structures were refined us- 
ing the FullProf Suite package 40 . 



III. RESULTS 

A. Nuclear structure of BiMnaMruO^ at room 
temperature 

We first consider the room temperature diffractogram 
of Fig. 2. The data analysis shows that the as- 
prepared samples are 96% pure or better with Mn3 04 
(haussmannite) as main impurity 41 . In agreement with 
Ref^, the room temperature structure was success- 
fully refined in the noncentro symmetric space group Im, 
with a=7.545(l) A, 6=7.362(1) A, c=7.536(l) A, (3 = 
91.178(2)°. In Fig. 2, we compare the observed and 
calculated diffractograms. The refined structural param- 
eters are reported in Table I, where we note a good agree- 
ment between observed and calculated intensities, with 
R p = 4.21%, wR p = 5.57%, and x 2 = 3.67. The salient 



features of the refined structure are seen by considering 
the simplified cubic 7ro3 structure. In this structure, the 
MnC>6 octahcdra are regular with the Mn 3+ ions (B sites) 
in the center of each octant of the unit cell; the square- 
coordinated Mn 3+ ions (A' sites) are located in the mid- 
dle of the cell axes; the Bi 3+ ions (^4 sites) occupy the cell 
corners and center; all twenty four oxygen atoms occupy 
the same site with high m.. symmetry. In the real Im 
structure of Table I, all sites split into sites with lower 
symmetry and lower multiplicity. For instance, the cen- 
trosymmetric eightfold B site of Mn with cubic .3. sym- 
metry splits into two noncentro symmetric fourfold sites 
Mn4 and Mn5 with translational .1. symmetry only. The 
displacements of the Mn 3+ and 2 ~ ions allowed by the 
Im symmetry lead to a macroscopic polarization of the 
crystal, P, approximately oriented along the ac-diagonal, 
in accord with the quasi isotropic properties of the ac- 
plane. On the other hand, the local dipoles (or micro- 
scopic polarization) of the MnO@ octahcdra have no con- 
straints because the symmetry of the Mn4 and Mn5 sites 
is purely translational. Thus, these dipoles have a non- 
zero component along the 6-axis. In Table II we report 
an estimate of the dipoles computed using a point-like 
charge model. 

In Table III we report the six Mn-0 bond lengths of 
the MnOg octahedra and note that they vary within a 
broad 1.840-2.295 A range. Hence, the description of the 



crystal field in terms of e s 



y 2 and 3z 2 — r 2 orbitals 



is no longer appropriate, which prevents a straightfor- 
ward application of the Goodenogh-Kanamori- Anderson 
(GKA) rules 4 ^ to estimate the sign and the strength of 
the exchange interaction between neighboring Mn 3+ ions. 
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FIG. 2. (Color online) Rietveld refinement of the neutron 
powder diffraction data of BiMn3Mn40i2 at room temper- 
ature. The wavelength used is A=1.225 A. Circles and con- 
tinuous line indicate the observed and calculated profiles, re- 
spectively. Vertical bars mark the positions of the calculated 
Bragg reflections. The continuous line at the bottom is the 
difference between observed and calculated patterns. The re- 
fined structural parameters are reported in Table I. 
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TABLE I. Refined structural parameters of BiMn3Mn40i2 at 300 K within the Im space group. Numbers in parentheses 
indicate statistical uncertainty. Atomic coordinates x, y and z are given in reduced lattice units. Statistical indicators of the 
refinement are reported at the bottom. 



Space Group: Im; q=7.545(l) A, b=7.362(l) A, c=7.536(l) A, g = 91.178(2)° 
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Reliability factors with all non-excluded points for pattern: 

_R-factors (not corrected for background): i? p =4.21, wR p =5.57, 7? e a:p=2.93, x 2= 3-62 
Conventional Rietveld ii-factors for pattern: i? p = 16.2, w_R p =16.2, R exp =8.50, x 2 — 3.62 
Global user-weigthed x 2 (Bragg contribution)=3.67 



a Anisotropic thermal factors (xlO 4 ): /8u=ll(8), /3 2 =215(25), /3 33 =53(16), 0i 2 =O, /3i 3 =10(8), /3 23 =0. 



TABLE II. Components along the a-, b- and c-axis of the local 
electric dipole d on the Mn06 octahedra estimated using a 
point-like charge model for the Mn 3+ and O 2- ions. Units 
are g e A, where q e denotes the electron charge. 
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TABLE III. Mn-0 distances in A for the octahedral Mn 3+ B 
sites. Numbers in parentheses indicate statistical uncertainty. 
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Another reason of inapplicability is the large buckling of 
the MnOg octahedra. Indeed, the Mn-O-Mn bond angle 
of ~ 135° falls in between the opposite limits of 90° and 
180° considered by the above rules. 



B. Low-temperature diffraction data: magnetic 
orderings of BiMn3Mn40i2 

We now consider the evolution of the diffraction data 
in the 1.5-100 K range (see Fig. 3). These data indi- 
cate that the symmetry of the nuclear structure is Im 
at all temperatures. At 55 and 25 K, two new sets of 
diffraction peaks show up, which are attributed to the 
two AFM orderings probed by previous magnetization 
measurements^. The indexing of the two sets is con- 
sistent with two distinct ki=(l/2,Q,l/2) and k 2 =(0,0,0) 
propagation vectors, respectively. This is similar to the 
case of the isovalent LaMii3Mn40i2, which exhibits an 
independent AFM ordering of the Mn 3+ B and A 1 sub- 
lattices at comparable temperatures Tn_b—78 K and 
r/v,A'=21 K, respectively. Since the exchange interac- 
tion between Mn 3+ ions is expected to be similar in 
BiMn3Mn40i2, for BiMn3Mri40i2 as well we propose 
the same scenario of two independent AFM orderings for 
the B and A' sublattices, with Tjv,b=55 K and TV,A' = 25 
K. The question of the absence of coupling between the 
two magnetic sublattices in quadruple perovskites has 
been discussed elsewher o 33 ! 38 . 

Fig. 3 shows that the AFM transition at 55 K is char- 
acterized by an unusual feature. The appearance of the 
magnetic peaks indexed by ki is followed at slightly lower 
temperature (50 K) by a sudden increase of intensity of 
certain Bragg peaks. Since the magnetic response of the 
system is AFM, this increase indicates that the AFM 
structure of the B sublattice contains a uniform com- 
ponent described by the k 2 propagation vector. On the 
other hand, the AFM ordering of the A' sublattice is 
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completely described by k2 and corresponds to the loss 
of the I centering (anti-body centering) of the lattice. 



C. Refinement of the magnetic structure of 
BiMn 3 Mn 4 Oi2 

To solve the magnetic structures of both sublattices, 
we determined the basis functions of the irreducible rep- 
resentations of the Im space group for the relevant prop- 
agation vectors and Wyckoff sites of the Mn 3+ ions (see 
Table IV). The T mag — 3Ti © decomposition into 
irreducible representations is obtained for the two Mn4 
and Mn5 sites of the B sublattice and for both ki and 
k2 vectors, whilst the T mag = Ti © 2T2 decomposition is 
obtained for the three Mnl, Mn2 and Mn3 sites of the A' 
sublattice and for the only relevant k2 vector. Thanks 
to the above analysis, we performed simulated annealing 
refinements 43 of the diffraction data with the following 
constraints: 1. the total magnetic moment has been set 
to zero for both A' and B sublattices, as the magneti- 
zation data of Ref. 2 — give evidence of an AFM response; 
2. we have assumed that the candidate structures follow 
either the Ti or T2 representation, not an admixture of 
the two. For the refinements, the two 1.5 and 35 K data 
set were used in order to solve separately the structures 



TABLE IV. Basis functions, xp, for the axial vectors associated 
with the irreducible representations Ti and F2 for the two 
different B and A 1 sites of the magnetic Mn 3+ ions and for 
the two different propagation vectors ki and k2 that describe 
the magnetic structure of BiMn3Mn40i2. 
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of the A' and B sublattices. Indeed, the ordering of the 
B sublattice is complete at 35 K, as apparent from the 
temperature dependence of the magnetic peak intensities 
in Fig. 3 and of the lattice parameters in Fig. 5. Also, 
35 K is sufficiently higher than the ordering temperature 
Tn,A' — 25 K of the A' sublattice, so the structure of the 
latter can be neglected. 
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FIG. 3. Temperature dependence of the powder neutron 
diffraction patterns of BiMn3Mn40i2- The wavelength used 
is A=2.43 A. Three magnetic transitions are visible at 55 (▼), 
50 (♦) and 25 K (V). The first two ones refer to the ordering 
of the magnetic _B-sublattice, whilst the third one refers to 
the ^'-sublattice. The (T) symbol labels the peaks associ- 
ated with the ki = (l/2, 0,1/2) propagation vector, whilst the 
(♦) symbol labels the Bragg peaks whose intensities exhibit a 
jump at 50 K and are therefore associated with the k2 = (0,0,0) 
propagation vector. The (*) symbol labels the peaks of the 
antiferromagnetic transition of the MnsCU impurity. Inset: 
detail of the behavior of the (101) and (HO)-(lOl) doublet 
Bragg peaks and of the (100) ± ki-(001) ± ki magnetic dou- 
blet peak. 



1. Magnetic structure of the B-sublattice 

As to the B sublattice, the magnetic structural models 
described by T\ yielded the best statistical results. The 
result of the refinement is shown in Fig. 4 and Table V, 
whilst the structure is illustrated in Fig. 5. We outline 
the characteristics of this structure by analyzing sepa- 
rately the two ki- and k2-components with moments par- 
allel and perpendicular to the ac-plane, respectively. The 
ki-component consists of an E-type AFM structure with 
moments in the ac plane Mb^c= 1-37 [Iq. This struc- 
ture can be viewed as zig-zag chains along the a- and 
c-axis with alternate antiparallel and parallel moments. 
We recall that LaMn3Mn40i2 displays a different C- 
type structure made of ferromagnetically (FM) coupled 
AFM ac-planes 34 . It is remarkable that two compounds 
with similar nuclear structures and isoelectronic proper- 
ties exhibit radically different magnetic orders. This is a 
first indication of the role played by the polar Bi 3+ ion 
on the magnetic ordering. As to the magnetic structure 
of the k2-component, Table V and Fig. 5 show that the 
moments are oriented along the 6-axis and that they are 
as large as Ms,h=l /is- In view of the discussion below, 
it is noted that the orientation depends on the crystal- 
lographic site, i.e. the moments of the Mn4 (Mn5) sites 
are parallel (antiparallel) to the 6-axis. In summary, the 
moments of the B sites result from an unusual super- 
position of two components Mg inc and Ms,6 with dis- 
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tinct ki and k2 propagation vectors. Hence, the magnetic 
structure can be viewed as a large and uniform modula- 
tion iMj,! of an _E-type structure. The total moment 
resulting from this superposition is Afe=l .70 hb, a much 
smaller value than the value of 4.0 hb expected for the 
high spin state (5=2) of the Mn 3+ ion. This discrep- 
ancy was reported for LaMn3Mn40i2 as well2i, thus 
confirming the inadequacy of a purely ionic model for 
the electronic states of manganese oxides with quadruple 
perovskite structure. 
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FIG. 4. (Color online) Rietveld refinement of the powder 
neutron diffraction data of the preceding figure. Top panel: 
refinement of the 35 K data. Circles and continuous line 
represent the observed and calculated diffraction patterns, 
respectively. Vertical bars mark the calculated position of 
the nuclear peaks (first row), the magnetic peaks of the ki- 
component of the Mn(B) structure appearing at 55 K (second 
row), the magnetic peaks of the k2-component of the Mn(5) 
structure appearing at 50 K (third row) and the nuclear peaks 
of the Mn304 impurity (fourth row). The continuous line at 
the bottom is the difference between observed and calculated 
patterns. Bottom panel: the same as for the top panel for 
the 1.5 K diffraction pattern. The additional row of ticks 
indicates the magnetic peaks of the Mn(A') structure. 



2. Magnetic structure of the A 1 -sublattice 

Also in the case of the AFM structure of the A' sublat- 
tice, the best structural model is described by the Ti rep- 
resentation. The magnetic moments and the structural 
model are reported in Table V and Fig. 6, respectively. 
The structure consists of an anti-body centered AFM or- 
dering with a moment Ma'—1-22 [Ib per site. The struc- 
ture turns out to be identical to that of NaMn3Mn40i2 33 
but different from that of LaMn3Mn40i2^- This con- 
firms the key role played by Bi 3+ in altering the stabil- 
ity of competing magnetic orderings in these two isoelec- 
tronic and almost isostructural compounds. Similarly to 
the previous case of the S-site, also the moment per A' 
site turns out to be much smaller than that expected 
according to a purely ionic model. 



IV. CORRELATION BETWEEN LATTICE 
POLARIZATION AND MAGNETIC ORDERING 

We should now provide a physical interpretation of the 
magnetic orderings observed. Our point concerns the 
symmetry properties of the Bi 3+ -induced lattice polariza- 
tion and of the AFM structure of the Mn 3+ B-sublattice 
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FIG. 5. (Color online) Schematic representation of the two 
components with propagation vectors ki = (l/2, 0,1/2) (left 
panel) and k2=(0,0,0) (right panel) which form the AFM 
structure of the Mn 3+ B-sublattice of BiMn 3 Mn 4 0i 2 . The 
purple balls indicate Mn 3+ ions. The purple arrows indicate 
the moments of the Mn 3+ ions, whilst the ± signs indicate 
the sign of the 6-component of the moments whilst the black 
(white) arrows (not in scale) indicate the b- (c-) components 
of the electric dipoles formed by the MnOe octahedra (the 
small a-component is not shown for simplicity). The mag- 
netic moments of the structure in the left panel are reversed 
in the adjacent nuclear cells along the a- and c-axis, as the 
structure is described by the ki propagation vector. This is 
not the case for the structure of the right panel, as this is 
described by k 2 =(0,0,0). 
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TABLE V. Magnetic structures of the B and A' sublattices of BiMnaMn40i2 at 1.5 K. The structures are obtained after 
refinement of the diffraction data of Figure 4. The nuclear unit cell parameters are: a= 7.537(1) A; 6=7.354(1) A; c=7.534(l) 
A; a, 7=90°; /3= 91.213(1)°. Magnetic moments are given in Cartesian coordinates and in Bohr magneton (ib units. Numbers 
in parentheses indicate statistical uncertainty. Approximate atom coordinates are given. The magnetic structure of both 
sublattices is described by the Fi representation, i.e. the moments (M x ,M y ,M z ) transform as ( — M x ,M y ,—M z ) upon application 
of the mirror operation (x,y,z) —¥ (x,—y,z) on atom coordinates (see text and Table IV). 





Mn 3+ B-sublattice (T N . B 


=55 K) ; propagation vector ki = (l/2 1/2) 




Atom coordinates 


M 


M x 


My 


M z 


Mn4 (1/4 3/4 1/4) 


1.37(1) 


-1.22(4) 





0.60(6) 


Mn4 (3/4 1/4 3/4) 


1.37(1) 


-1.22(4) 





0.60(6) 


Mn5 (1/4 3/4 3/4) 


1.37(1) 


1.22(4) 





-0.60(6) 


Mn5 (3/4 1/4 1/4) 


1.37(1) 


-1.22(4) 





0.60(6) 


Mn 3+ B-sublattice (Tjv,s=50 K); propagation vector k 2 = (0 0) 


Mn4 (1/4 3/4 1/4) 


0.99(2) 





0.99(2) 





Mn4 (3/4 1/4 3/4) 


0.99(2) 





0.99(2) 





Mn5 (1/4 3/4 3/4) 


0.99(2) 





-0.99(2) 





Mn5 (3/4 1/4 1/4) 


0.99(2) 





-0.99(2) 





Mn a+ ^'-sublattice (Tjv iA /=25 K); propagation vector k 2 = (0 0) 


Mnl (0 1/2) 


1.22(1) 





1.22(1) 





Mnl (1/2 1/2 0) 


1.22(1) 





-1.22(1) 





Mn2 (1/2 1/2) 


1.22(1) 





-1.22(1) 





Mn2 (0 1/2 0) 


1.22(1) 





1.22(1) 





Mn3 (1/2 0) 


1.22(1) 





1.22(1) 





Mn3 (0 1/2 1/2) 


1.22(1) 





-1.22(1) 






reported in Fig. 5 and in Table V. We show below that 
these properties point to a large magnetoelectric coupling 
in BiMn3Mn40i2. We first note that the order param- 
eters of the polarization and of the 6-component, Mb, 
of the magnetic moments M are both described by the 
k2 propagation vector. This suggests that Mb is a uni- 
form modulation of the E-type structure induced by the 
noncentrosymmetric displacements of the Mn 3+ and 2 ~ 
ions via a magnetoelastic coupling. To show the validity 
of this scenario, we recall that the dominant contribution 
of this coupling to the Landau free energy is of the form 




FIG. 6. The same as in the preceding Figure for the Mn 3+ 
A'-sublattice. 



/ = —^^CijMkMi, where j3 is a fourth-rank coupling 
tensor, e is the strain tensor and i,j,k,l are the indices 
of the a, b and c-axis^ 4 -. In our case, the coupling arises 
from the internal strain associated with the dipoles of 
the Mn0 6 octahedra (see Table II and Fig. 5). Indeed, 
this field uniformly alters the electrostatic potential and 
the spin-orbit interaction of the Mn 3+ ion as well as the 
exchange integrals of the Mn-O-Mn bonds. This prob- 
lem was extensively treated by Callen and Calient, who 
considered both the external strain and the internal one 
caused by the atomic displacements. We recall that, by 
symmetry, the spins can couple only to strain modes sym- 
metric under spatial inversion. In our case, due to the 
mirror ac plane of our pseudocubic Im symmetry, the 
only relevant mode is the Ti mode symmetric under the 
y — )• — y inversion. This observation confirms the validity 
of our scenario of magnetoelastic coupling between Mb 
and the polarization-induced strain, since the Ti repre- 
sentation also describes Mb- The reason is that the free 
energy contains the fully symmetric representation to be 
obtained from the direct product of the two representa- 
tions of the strain and spin basis function o 44 ' 45 . Thus, the 
two representations must be equivalent, as in our case. 
The existence of a coupling is confirmed by the fact that 
an internal strain along the 6-axis does indeed exist, as 
the dipoles of the MnOe octahedra have a sizable compo- 
nent along the 6-axis. Owing to the mirror ac-plane, the 
atomic displacements associated with the dipoles create 
a modulation of the Mn-Mn distance along the 6-axis (i.e. 
an internal strain ebb), the short (long) distance being as 
small (large) as 3.490 (3.858) A. 

The above considerations also account qualitatively for 
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the stability of the M& structure. Indeed, the moments 
are bound to be parallel along the 6-axis because of the r x 
symmetry of the strain field. In addition, the moments 
on the Mn4 and Mn5 sites tend to be antiparallel as the 
exchange interaction is predominantly AFM within the 
ac-plane 3 -^. Both characteristics are indeed found in our 
experiment. A quantitative account of the stability of 
the Mb structure goes beyond the scope of this work; it 
may be provided by ah initio calculations in the local 
spin density approximation that include the on-site re- 
pulsion between d-electrons or ad hoc Hamiltonians that 
include the charge, orbital and spin degrees of freedom. 
These methods have been applied to LaMnsMruiOia — , 
i?Mn03^2r— and -RM^Os^S among others, and an ap- 
plication may be envisaged in our case as well. 

Our picture of magnetoelastic coupling along the 6-axis 
also accounts for the temperature-dependence of the unit 
cell parameters reported in Fig. 7. One notes that the 
AFM ordering of the B sites at 55 K is accompanied by a 
sizable expansion of both a- and c-axes, whilst the 6-axis 
displays no anomaly. This indicates a magnetostriction 




T(K) 



FIG. 7. Temperature dependence of the lattice parameters of 
BiMn3Mn40i2. The data are obtained from the refinement of 
the powder neutron diffractograms of Fig. 2. Note the large 
magnetostriction along the a- and c-axis at the antiferromag- 
netic ordering at 55 K of the Mn 3+ ions in the octahedral 
B-sublattice (see text). 



only in the ac-plane, which is explained as follows. The 
variations of the unit cell parameters reflect the external 
strain induced by the macroscopic polarization, P. As 
mentioned before, in BiMn3Mn40i2, P is approximately 
oriented along the ac diagonal, while its 6-component is 
identically zero. Therefore, we expect a magnetostric- 
tion only along this diagonal, in agreement with Fig. 
7. Similar arguments account for the absence of mag- 
netostriction at the AFM transition of the Mn 3+ ions at 
the A' sites at T W .A'=25 K. Indeed, the ^-component of 
the dipole at this site is identically zero because of the 
m symmetry of these sites. Thus, the only symmetric 
internal strain mode Ti available for the magnetoelastic 
coupling vanishes in this case. 



V. CONCLUSIONS 



In conclusion, we found that the magnetic order of the 
Mn 3+ ions in the quadruple perovskite BiMn3Mn40i2 is 
strongly modulated by the internal strain field induced 
by the polar Bi 3+ ion. This observation suggests a mech- 
anism of magnetoelectric coupling between Bi 3+ and 
Mn 3+ ions driven by the magnetoelastic coupling due 
to the internal strain field, in agreement with a previous 
observation of giant magnetoelastic effect in the multifer- 
roic i?Mn03 system^. Our analysis indicates that this 
mechanism is effective in BiMn3Mn40i2 thanks to two 
features of quadruple perovskites: a) the high symme- 
try of the oxygen atoms surrounding the Bi 3+ ion, which 
limits structural distortions, thus preventing the release 
of the strain; b) the absence of oxygen disorder. Further 
experimental and theoretical work is needed to test the 
tunability of the magnetic order of BiMii3Mn40i2 by an 
electric field, a prerequisite for applications. 

In the light of the above considerations, our result 
demonstrates that a large magnetoelectric coupling be- 
tween magnetic and polar ions can be achieved in proper 
ferroelectrics. The coupling mechanism which accounts 
for our result suggests that site symmetry is an effective 
control parameter of the internal strain responsible for 
the coupling. This conclusion suggests a novel concept 
of internal strain engineering for the effective design of 
multifcrroics. 
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